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ABSTRACT 

We study the gravitational Faraday rotation, on linearly polarized light rays emitted 
by a pulsar, orbiting another compact object. We relate the rotation angle to the 
orbital phase of the emitting pulsar, as well as to other parameters describing its 
orbit and the orientation of the angular momentum of the binary companion. We 
give numerical estimates of the effect for the double-pulsar system PSR J0737-3039, 
and we note that the expected magnitude is exceedingly small, making the effect 
unlikely to be observed with present technology. It is however interesting per se, since 
in this phenomenon, gravito-magnetism plays a leading role, unlike what happens, for 
instance, when studying light bending or gravitational time delay, where it appears as 
a correction to the gravito-electric contribution. Also, we envisage the possibility that 
this effect could be relevant, at least in principle, for a pulsar orbiting a non charged 
black-hole. 

Key words: gravitation, relativity, polarization, pulsars: general, pulsars: individual 
(PSR J0737-3039) 



1 INTRODUCTION 

The physical properties of the rece ntly discovered double - 
pulsar sys t em P SR J0737-3039 (iBurgav et al.l J2003D . 
iLvne et alJ i2004) ') make this system a rare laboratory for 
testing relativistic gravity. Namely, the system is composed 
of a 22-ms pulsar, PSR J0737-3039 A, in a slightly eccentric 
(e ~ 0.088) 2.4/1 orbit, together with PSR J0737-3039 B, 
whose period is 2.7 s: thanks to the high mean orbital speed 
{v ~ 0.001c) this system is highly relativistic, and allows the 
detection of four post-Keplerian parameters, thus accurately 
constraining the masses of the two stars and providing also 
string ent tests of General Relativity (GR) ijKramer et alJ 

Actually, excluding cosmological issues, GR has passed 
all observational tests with excellent results ( Will (2006)). 
However, we must remember that most of the tests of GR 
come from Solar System experiments, where the gravita- 
tional field is in the "weak" regime. On the other hand, it is 
expected that deviations from GR can occur in the "strong" 
field regime: hence, the Solar System experiments are inade- 
quate to this end. On the contrary, the strong gravitational 
field is best tested by means of pulsars. Among the post- 
Newto nian effects that should be relevant in strong field 
regime l|WilJ (|l993f li. there are the so called gravito -magnetic 
(GM) effects, which arise from the rotation of the sources 
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of the gravitational field: these effects are generally small, 
since the gravitational coupling with the angular momen- 
tum of the source is much weaker than the coupling with 
mass alone (the so called gravito- electric interaction), and 
this m a kes their detection very diffic ult ('see lMashhoon et alJ 
l^2001^ : lRuggiero fc Tarta dia' ('20Qj)). 

In previous pa pers (Tart aglia fc Ruggierd i2005l) : 
iRuggiero fc Tartaglial ifioOiifl i. we pointed out that binary 
pulsar systems in general and, in particular, the double- 
pulsar PSR J0737-3039, could be very interesting for test- 
ing the GM effects. The relevance of GM effects in bi- 
nary pulsar systems was investigated in the past and also 
more recently by many authors, who also pointed out 
the difficulties that the detection of these effects must 
face I Doroshen ko fc K opeikinI lll995l): iLaguna fc Wolszczari 
1*1997): IWex fc Ko peikin ('1999'): iKopeikin fc Mashhoon 
(200l: lRafikov fc L ai ( 2006a b)), mainly due to their smaU- 
ness, which is overwhelmed by other dominant effects. 

Here, we would like to study a physical phenomenon 
where the gravito-magnetic interaction has a leading role, 
and which can be tested, at least in principle, in binary pul- 
sar systems, i.e. systems where a pulsar is orbiting a com- 
pact object, such as another pulsar, as in PSR J0737-3039, 
or, hopefully, a black hole. The phenomenon we are going 
to deal with is the gravitational Faraday rotation. 

The polarization vector of an electromagnetic (e.m.) 
wave propagating in a gravitational field undergoes a change 
of direction, as a result of the deflection of e.m. signals, and, 
in addition a rotation around the propagation vector. As 
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for the latter, in particular, it has been studied and dis- 
cussed by many aut hors, in the context of GR (see, for early 
works the paper by IPIe bans kj fW6&) and, for more recent 
references the paper bv iSereng (.2005..) ). Actually this gravi- 
tational rotation of the plane of polarization, for linearly po- 
larized light rays, is the gravitational analogue of the electro- 
magnetic Faraday effect, that is the rotation which a light 
ray undergoes when it propagates in a plasma under the 
action of a magnetic field, the role of the magnetic field be- 
ing played, in this case, by the gravito-magnetic field, ft is 
possible to show that, for astrophysical interesting cases, a 
spherical and non rotating source of gravitational field does 
not cause a net rotation of the polarizat i on vec tor jPlebanskil 
Jl960h : lDver fc Shaveil Jl992l) : lFaraonil lll993l) ). while things 
change when rotating sources come into play. Indeed, the 
gravitational Faraday rotation is a purely gravito-magnetic 
phenomenon (Sereno ( 2005)). Actually, also the case of rota- 
tion is somewhat su btle, since in the l owest order of the weak 
field approximation. IPlebanskl (Il960l) concluded that the ro- 
tation of the polarization plane occurs only when the light 
rays penetrate into a rotating distribution of matter. More 
accurate calculations, that take into account the deflection 
of the light rays and higher-order terms in the weak field 
approximation, show that rotation o ccurs outside t h e hor i- 
zon of a rotating (Kerr) black ho le lllshihara et alJ il98d) : 
lNouri-Zono3 Jl999l) : ISerend ^2005^ ). 



The phenomena we try and study here could be partic- 
ularly relevant in the case of X-ray emission from accretion 
disks or other thick material envelopes surrounding rotat- 
ing magnetized compact objects and specifically black holes. 
Indeed the intrinsic birefringence of the plasma in the mag- 
netic field, together with the general relativistic effects due 
to the curvature in proximity of horizons, can be the pri- 
mary source of polarization and of the rotation of the po- 
larization plane of th e emitted short wavelength rad iation 
JConnors et alJ Jl980l) : lBroderick fc BlandfordH2004) ). 



Under the hypothesis that, to a reasonable order of ap- 
proximation, the space-time around a rotating compact ob- 
ject (i.e. a pulsar, or a black hole) can be described by the 
Kerr metric, here we investigate the gravitational Faraday 
rotation induced on the light rays emitted by a pulsar or- 
biting around it. We study the effect and give estimates of 
its magnitude, and evaluate the possibility of detecting it. 



The paper is organized as follows: after briefly reviewing 
the basics of propagation of e.m. signals in a gravitational 
field that lead to the gravitational Faraday rotation in Sec. 
|5| we apply the formalism to the case of the signals emitted 
by a pulsar orbiting a compact object in Sec.|21 The estimate 
of the magnitude of the effect, together with the evaluation 
of the possibility of detecting it, is discussed in Sec.|ll while 
the conclusions are summarized in Sec.|3| 



2 PROPAGATION OF E.M. SIGNALS IN A 
GRAVITATIONAL FIELD AND THE 
GRAVITATIONAL FARADAY EFFECT 

The propagation of e.m. signals in a gravitational field, can 
be studied by writi ng the Maxwell equat ions in the Lorentz 
gauge; on following iFavos fc Llosal I^S^), we may write ^: 

V,,V"A'' = 0, (f) 



Vo,A" = 0. (2) 

This is done considering two basic approximations; (i) 
the electromagnetic field is weak enough not to affect the 
gravitational field; (ii) the geometric optics limit is adopted, 
where it is assumed that the wavelength A of the e.m. waves 
is much smaller than L, where L — min{jC, R), £, being 
a typical distance where amplitude, polarization and wave- 
length vary appreciably and R is the order of magnitude of 
the (spatial) curvature radius. Then, we may introduce the 
wave four- vector k" and polarization four- vector f", which 
fulfill the following relations 

k"k^ = 0, fc"/c = 0, ,ru = 1- (3) 

Light rays follow null geodesies, whose tangent vector 
is k"; furthermore, the polarization four vector undergoes 
parallel transport along these geodesies. Consequently, we 
may write 

rVcfc*^ = 0, (4) 

k"Vc,f = 0. (5) 

The changes induced by the gravitational fleld on the 
wave and polarization four-vectors, are then obtained inte- 
grating the equations (I4I5II . Actually, this problem can be 
suitably studied in the f-|-3 projection formalism, applied to 
the stationary space-time describing the gravitational fleld 
in which w aves propagate. This approach is outlined by 
iFavos fc L losa (f982); Nouri-Zonoz (i999): here, we briefly 
recall its foundations. 

Let M be the manifold which is the model of our space- 
time, endowed with the metric gap, together with a metric 
connection F. Since the space-time is stationary, there exists 
a time- like Killing vector field ^j. Then, one can show that, 
for each p £ A4, there is a three-dimensional manifold E3, 
defined by the smooth map 

* : ^ E3, (6) 

^ — *I'(p) being the orbit of the Killing field passing 
through p. We can introduce a coordinate system adapted 
to the congruence = dt, and let ytj be the projected 
three-dimensional metric of E3. Suitable differential oper- 
ators can be also defined on E3, starting from jij, in such 

^ Greek indices run from to 3, Latin indices run from 1 to 3; 
the space-time metric has signature (1, — 1, — 1, — 1), and we use 
units such that G=c=l; boldface letters refer to three dimensional 
vectors. 
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a way that a covariant derivative V for th e three vectors 
belong i nR to S3 is defined (see, for details [ F avos fc Llosal 
l)l982l) : iNouri-ZonoJ l)l99d ): iLandau fc LifshitJ l|l97ll) and 

references therein). 

If we write the metric of a stationary space-time in the 
form 



where 



and 



= , h, = goo, 

goo 



(7) 
(8) 



dr 



-gii + 



500 



dx^ dx-' = 'yijdx^ dx-' , (9) 



it is possible to introduce the gravito-electric and gravito- 
magnetic fields, given respectively by 

E, = -Vh}'^, (10) 



Bg^V X A. 



(11) 



In this context, on the basis of an orthogonal decom- 
position in adapted coordinates, the three- vectors A; and /, 
representing the projection of k" , on E3, can be taken 
to be equivalent to the contravariant components of k" and 

respectively. Then, from equations 14151 . it is possible 
to write the three-dimensional equations 



Vfefc = L X fe + [Eg ■ k) k 



(Eg,Bg,k,f) (see iFavo 



fc Llosal 



(12) 

(13) 
1982 



where L = L 
iNouri-Zon ol fm\ 

Eqs. I12I13II show, in particular, that the polarization 
vector rotates with angular velocity L along the projected 
geodesic. By carrying out explicit calculations, the angle of 
rotation around fc which expresses the gravitational Faraday 
rotation is given by 



eta = -- 



VhBg ■ dl, 



(14) 



where the integral is evaluated from the source (emission 
point) to the observer, along the null geodesies whose tan- 
gent (spatial) vector is fc, such that dl — kdl. 

It is useful to recall he re the expression of the el ectro- 
magnetic Faraday rotation ijLorimer fc Krameil (|2(i 



O^Faraday — 



nB\\ dL 



(15) 



where we have the physical quantities that describe the prop- 
erties of the plasma where the wave propagates (n is the 
refraction index), the angular frequency ui of the radiation 
and -By , which is the component of the magnetic field along 
the line of sight. On comparing (1141 and (1151 we see that 
both depend on the component of the gravito-magnetic or 
magnetic field along the line of sight. There is a formal anal- 
ogy, even though the effects are physically different since the 



gravitational Faraday rotation is a purely geometric effect, 
while the traditional Faraday rotation depends on the fre- 
quency of the light rays. 

Formula (1141 shows also that the gravitational Faraday 
rotation is an effect due to mass-current, i.e. to the rota- 
tion of the source of the gravitational field, so it is a purely 
gravito-magnetic phenomenon, and it is null in space-times 
around non rotating sources. 

As we said before, this effect was investigated in the 
past by many authors, in suitable weak field approximations. 
We are interested in the effect on light rays emitted by a 
pulsar orbiting a compact o bject: henc e , we c annot apply 
here calculations outlined in (piebanskil lll960f) , since their 
level of approximation is too low, and gives no rotation of the 
polarization vector unless light rays penetrate into rotating 
distributions of matter. 

A better level o f approximation wa s considered b y 
llshihara et ai] Jl98d) : iNouri-Zonoj |l99i); ISerend tOO^ . 
who studied propagation of light rays outside the horizon of 
a rotating Kerr black hole. These authors agree on the right 
order of approximation needed, however there is some dis- 
agreement on numerical values. In particular, iercno (200|) 
obtained for the gravitational rotation angle the expression 



TVMJ« 



(16) 



where /i is a parameter that quantifies the contribution of 
angular momentum to space-time curvature (^ = 1 in GR), 
M is the mass of the source of gravitational field and J\\ is 
the component of its angular momentum along the line of 
sigh, ^ is the projection in the plane of sight of the distance 
between the emission point and the source of the gravita- 
tional field. 

In the following Section, under the hypothesis that, to 
a reasonable order of approximation, the space-time around 
a rotating compact object (i.e. a pulsar, or a black hole) can 
be described by the Kerr metric, we apply (1161 in order to 
study the gravitational Faraday rotation on the light rays 
emitted by a pulsar orbiting the compact object. 

Before going on, we notice that a better (or more re- 
alistic) approximation for describing the space-time around 
a rotating compact object such as a pulsar, is given by the 
metric obtained bv iHartle fc Thornd il968ll . Since a neu- 
tron star's moment of inertia is larger than the one of a 
black hole of similar mass, the Hartle-Thorne metric could 
lead to a slightly larger gravitational Faraday rotation, how- 
ever we do not expect that this can significantly change the 
order of magnitude of the effect. 



3 GRAVITATIONAL FARADAY EFFECT IN A 
BINARY PULSAR SYSTEM 

We use the notation of Figure Q for the description of the 
pulsar orbi t around the compact binary companion (see 
IStraumannI (120041) '). We choose a first set of Cartesian co- 
ordinates {x,y,z}, with origin in the center of mass of the 
binary system, and such that the line of sight is parallel 
to the 2 axis. Then, we introduce another set of Cartesian 
coordinates {X,Y,Z}, with the same origin: the X axis is 
directed along the ascending node, the Z axis is perpendic- 
ular to the orbital plane. The angle between the x and X 
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Figure 2. Notation and conventions for describing the pulsar or- 
bit and the spin axis of the binary companion. The spin axis J ori- 
entation is determined by the angles Xc,t]c- The system {x,y,x} 
is such that xy is the plane of sight, the {X, Y, Z} system is such 
that XY is the orbital plane. 



Figure 1. Notation used for describing the pulsar orbit. 

axes is Q, the longitude of the ascending node, while the 
angle between the z and Z axes is i, the inclination of the 
orbital plane. 

Let be the vector describing the orbit of 

the pulsar (here and henceforth, the suffix p refers to the 
emitting pulsar, while c refers to its binary companion): it 
is described by 

Xp — cos(tj + ^p)X + sin(cij + ip)Y , (17) 

in terms of the argument of the periastron, ui, and the true 
anomaly, ip. Let us pose 

e^w + ip, (18) 

for the sake of simplicity. Then, we use the notation 

r = Xp~Xc (19) 

to describe the position of the pulsar with respect to its com- 
panion, and we remember that we have, for the Keplerian 
problem 



a(l 



1 + e cos (f 



(20) 



where a is the semi-major axis of the relative motion and e 
is the eccentricity. The astronomical elements Q, i, ui, a, e 
represent the Keplerian parameters. 

Then, let J be the angular momentum of the binary 
companion: its direction is determined by the angle Ac == 
z ■ J, and rjc between the projection of the J onto the plane 
of the sky and the ascending node of the binary orbit (see 
Figure I^J . 

As a consequence, we may write 

J — sin Ac cos {^1 + rjc) X + sin Ac sin (fi + 77c) y + cos AcZ. 

(21) 



The unit vector along the line of sight is n = — z, so 
that we have 

J|l =|J|(-cosAc). (22) 

Furthermore since ^ — \n x {r x n) \, from I|17I19|I and 

n = — cos iZ — sin iY , (23) 

we have 

■ (24) 



^ = r\/ cos^ 9 + sin^ 9 cos^ 



Consequently, the rotation angle of the polarization vec- 
tor llfcil l. setting /i = 1 as in GR, becomes 



TT M\J\ 



cos Ac 



4 (cos2 6»-^sin2 6'cos2i) 



3/2 ■ 



(25) 



or, explicitly writing the equation of the orbit 12011 and re- 
introducing the mean anomaly 



TvG^ M\J\ 



cos Ac (1 + e cos ip) 



4 c5 a3(l - e2)3 (cos2 + sin^ {u + <p) cos^ 



,3/2 ' 



(26) 

in physical units. 

We see that the rotation angle depends on the orbital 
phase: this fact, could lead, in principle, to a way of detecting 
it. 



4 DISCUSSION 

In order to quantitatively evaluate the magnitude of the 
gravitational Faraday rotation for light rays emitted by a 
pulsar in a binary system, it is useful to write (I26II in the 
form 



(27) 
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Figure 3. Variation of the phase F{ip), along the pulsar orbit: 
above, the inclination of the orbit is the one measured in the 
double-pulsar PSR J0737-3039; below, two smaller values have 
been chosen. 



where 



A 



l9l 

4 c5 



M\J\ cos Ac 
a3(l - e2)3 



(28) 



might be thought of as an "amphtude", independent of the 
orbital phase, while 



(1 + ecos ip) 



V 3/2 



(29) 



(cos2 [uj + (f) + sin^ (lu + (p) cos^ i)" 

is a "phase" factor, explicitly depending on the anomaly ip. 

For estimating A, we consider, for i nstanc e, data coming 
from PSR J0737-3039 iKramer et al.l ll2005ll '): accordingly, 
a = 9 X 10^ Km,e = 0.088, = 75.85°, sin i = 0.995, M = 
1.25OM0 (we considers light rays emitted by PSR J0737- 
3039 A, while it orbits PSR J0737-3039 B), and estimate 
I J| ~ \Jq\. Since we don't know the orientation of the spin 
of PSR J0737-3039 B, we may suppose, for the sake of sim- 
plicity cos Ac ~ 1, which is the most favorable situation. 
Consequently, we obtain 

-18. 



A ~ 0.94 X 10"'Vad. 



(30) 



The function F{ifi) is plotted in Figure |3 for different 
values of the inclination. We see that, for orbits with great 
inclination, the effect grows as far as cu + ip ^ 7i"/2, Stt/2, i.e. 
close to conjunction. 



We may say that, provided that Jy 7^ 0, the effect is big- 
ger for small orbits with great eccentricity and inclination. 
Actually, even though the peculiar dependence on the or- 
bital phase may give a chance to measure the rotation angle, 
in very favorable conditions, the overall effect remains very 
small, and furthermore, it is overwhelmed by the electromag- 
netic Faraday rotation (see, for instance, iDemorest et alJ 
i2004f) . and references therein). 

However, it is interesting to point out that since the 
electromagnetic Faraday rotation 1151 1 is proportional to the 
square of the wavelength of the radiation (inversely propor- 
tional to LJ^), the situation could be different in the case of 
short wavelength emission, such as X-rays. 

To give a quantitative evaluation, we recall that the 
pol arization position angle (PP A) can be written in the form 
(see iLorimer fc Krameil ll2005h ) 



aFPA(A) 



^Faraday 



y X RM, 



(31) 



where A is the wavelength, and the rotation measure RM is 
defined by 



RM 



nB\\dl. 



(32) 



The rotation measures of the two pulsars A and B in the 
double s ystem PSR J0737-303 9 have been recently mea- 
sured bv lPemorest et all i2004l) . and the corresponding val- 
ues are —112.3 ± 1.5 rad and —118 ± 12 rad , 
respectively. On using these values, we may estimate the 
PPA of the X-ray emission from PSR J0737-3039 A. On 
the ground of the recent observations of the X-ray emis- 
sions f rom the double pulsar perfor med bv lMcLauehlin et alJ 
i2004f) and lCampana et alJ i2004f) . we consider emission at 
1 keV{\ = 1.24 X 10"^m), which gives 



appA{\ = 1.24 X IQ'^m) ~ 1.72 x 10"^>ad. 



(33) 



In other words, if the other physical conditions are the 
same, the electromagnetic Faraday effect for X-rays is some 
17 orders of magnitude smaller than for radio wavelengths: 
in these conditions, the purely geometrical frequency inde- 
pendent rotation, i.e. the gravitational Faraday rotation, can 
become relevant. 



5 CONCLUSIONS 

We have studied the gravitational Faraday rotation on hght 
rays emitted by a pulsar, orbiting another compact object, 
such as a pulsar or a black-hole, whose gravitational field, 
to a reasonable order of approximation, can be described by 
the Kerr metric, of which we have considered the suitable 
weak-field limit. 

In this phenomenon, gravito-magnetism plays a leading 
role, unlike what happens, for instance, when studying hght 
bending or gravitational time delay, where it appears as a 
correction to the gravito-electric contribution. 

The effect is a purely geometric one, contrary to the 
electromagnetic Faraday rotation, which depends on the fre- 
quency of the hght rays; furthermore, it depends on the com- 
ponent of the angular momentum of the source of the grav- 
itational field along the line of sight. So, when the latter is 
not null, a net rotation of the polarization plane is expected. 

We have obtained a formula which relates the rotation 
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angle to the orbital phase of the emitting pulsar, as well as 
to the parameters describing its orbit and to the orientation 
of the angular momentum of the binary companion. 

The numerical estimates of the magnitude of this ef- 
fect, in the context of the double-pulsar system PSR J0737- 
3039, show that at radio wavelengths it is overwhelmingly 
dominated by the electromagnetic Faraday rotation. Never- 
theless, as a result of it's achromatic nature, gravitational 
Faraday rotation will dominate electromagnetic Faraday ro- 
tation at X-ray energies. However, in practice, due too its 
exceedingly small nature, it is unlikely to be observed in 
pulsar systems. 

Anyway, this effect could be of some relevance 
for pulsar-blac k hole binaries, or accreting black holes 
iConnors et alJ dlOSCf) '). In such systems, peculiar modula- 
tion of polarization orientation during the orbital motion 
may be evidence that an unseen companion is a black hole. 
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